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Abstract 



The nonlocal gauge invariant mass operator Tr J d'^x Fp,^{D^)~^ F^i,i, is investigated in Yang-Mills theories in 
the maximal Abelian gauge. By means of the introduction of auxiliary fields a local action is achieved, enabling 
us to use the algebraic renormalization in order to prove the renormalizability of the resulting local model to 
all orders of perturbation theory. 

I> 

d 1 Introduction 

^ ■ 

o : 

■ One of the major open problems in quantum field theory is the understanding of nonabelian gauge theories, and 
^ , consequently of quantum chromodynamics (QCD), in the infrared regime. The confinement phenomenon of quarks 
and gluons is not yet clearly established from the theoretical point of view and still waits for a satisfactory expla- 
nation. 



The Yang-Mills (YM) theories are described by the following Euclidean action 

5ym = jyAF^tF^t, (1) 

where is the field strength 

< = d,At - d^Al -f .g/-^^^A^ A? . (2) 

Here f^^'-^ are the structure constants of the gauge group SU{N) with A = 1, . . . ,7V^ — 1, and g is the coupling 
constant. At high energies, the running coupling constant is sufficiently small to allow for a perturbative descrip- 
tion, as expressed by the asymptotic freedom [HE]- However, when one lowers the energy, the running coupling 
constant grows, causing perturbation theory to fail, so that nonperturbative techniques are required. 
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To deal with this problem, different approaches have been considered. For example, in the Landau gauge, several 
analytic results have been obtained on the infrared behavior of the propagators of the theory, see for instance 
[SI m [5l [6l [71 m m Uni HT]. in this gauge, lattice simulations have confirmed an infrared suppressed gluon prop- 
agator exhibiting violation of positivity [HI [T31 dH [121 [TB], a feature interpreted as a signal of confinement. In 
particular, the fitting of the lattice data for the gluon propagator are usually accomplished with the aid of several 
mass parameters [T^ [T71 [TS], whose theoretical interpretation is still under investigation. So far, two possible 
origins could be suggested for such massive parameters, namely: the existence of the Gribov copies [3] and the 
condensation of suitable dimension two operators built up with gluon and ghost fields [191 (IHl [H] • In this work we 
shall focus on dimension two operators built up with gluon fields only, see [22l [23] for a recent discussion of ghost 
condensation in the Landau and maximal Abelian gauges. 

The massive Gribov parameter 7, which is fixed by a gap equation, follows from the restriction of the domain 
of integration in the Feynman path integral to the so called Gribov region [31 [31 [S] . This restriction is needed in 
order to deal with the problem of the Gribov copies and may be implemented through the introduction in the 
YM action of a nonlocal term, known as the Zwanziger horizon function ^4]. However, there is still room for 
additional mass parameters. Therefore, the possibility of the condensation of dimension two operators, giving rise 
to a dynamically generated mass for the gluons, has been taken into account. Let us also mention that, besides 
the lattice data, the introduction of an effective gluon mass turns out to be useful also from the phenomenological 
point of view, see for example [24l [25l [26] . In particular, in the Landau gauge, the dimension two operator A^A^ 
was proven to be multiplicatively renormalizable to all orders [23 . In [21] [28l |29l [30] , an effective potential for 
A^A^ was constructed. The formation of a nonvanishing condensate {A^A^)^ resulting in a dynamical effective 
gluon mass, turned out to be energetically favored. 

Besides the Landau gauge, other gauge fixings have been considered. We mention the recent analysis of dimension 
two operators in the Curd- Ferrari [31] and general linear covariant gauges [32l [33] . In [31] , an effective potential 
was constructed for the on-shell BRST invariant operator [jA^A^ + ac^c"^) in the Curci-Ferrari gauge while, in 
[32] [33], a detailed study of the already mentioned operator A^A"^ was performed in the linear covariant gauges. 
As a result, in [33], it was shown that the gluons do acquire a dynamical mass since the formation of (A^A^j^) 
is energetically favored. Another interesting gauge which is receiving increasing attention in the last years is the 
maximal Abelian gauge (MAG) [331 [351 [3B] . Several results have already been estabhshed in this gauge, both from 
theoretical [371 [351 [Ml [301 [U and lattice [31 [31 [33] points of view. The MAG is well suited for the study of special 
aspects of infrared QCD and color confinement as, for instance, the dual superconductivity and the so called Abelian 
dominance. The dual superconductivity mechanism [35] [3S1 [37] asserts that the low energy regime of Yang-Mills 
theories should exhibit monopoles as vacuum configurations. The condensation of these magnetic charges might 
give rise to a dual Meissner effect in the chromoelectric sector. As for the Abelian dominance hypothesis [3^, the 
infrared limit of QCD should be described by an effective theory constructed only from Abelian degrees of freedom, 
identified with the diagonal components of the gauge field, corresponding to the generators of the Cartan subgroup 
of the gauge group. Lattice numerical simulations in the MAG have reported significant differences between the 
diagonal and and off-diagonal components of the gluon propagator [311 [331 [33]- In particular, the off-diagonal 
gluon propagator displays a mass greater than that reported for the diagonal component, corroborating in fact the 
Abelian dominance hypothesis. In an attempt to understand those lattice results in the MAG, also here the con- 
densation of dimension two operators has been considered. In [371 149j . a dynamical mass generation mechanism for 
the off-diagonal gluons was proposed in the MAG, by means of the condensation of the operator {^A'^A^ + ac°c°)0. 

As the reader may have noticed, the dimension two operators mentioned above are gauge dependent, being related 
to specific choices of the gauge fixing. This is a consequence of the fact that a local gauge invariant dimension two 
operator is not available in YM theories. Still, the condensation of these operators might be taken as evidence 
in favor of the existence of a more fundamental gauge invariant operator. However, willing to preserve gauge 
invariance, we are led to give up of the locality requirement. The price one has to pay for that is that nonlocal 
operators are difhcult to be handled within a consistent renormalizable framework. So far, several possibilities have 
been considered. The first proposal for a condensate of dimension two was made by |19[ I20j. who considered the 

^Here the index a runs only on the off-diagonal components, see the beginning of the next section for the notations. 



2 



nonlocal gauge invariant operator A^^^, obtained by minimizing the operator A^A^ along the gauge orbit, namely 



^min = min 



jd\Tr{A^^)\ (3) 



where U represents an element of the gauge group SU{N), A^ = UA^U~^ + iUd^j^U^^. However, for a generic 
choice of the gauge fixing, the operator A^i„ proves to be very difficult to be handled at quantum level. Expanding 
A^jjj in a power series in the gauge field, see for example [50j . one obtains 



^min I d X 



^J^\ aA „fABC ((^c^aa\ ( 'Lf>aB 



0{A^). (4) 



The series ([4]) consists of an infinite number of nonlocal terms. So far, a consistent treatment for A^jj^ has been 
achieved only in the Landau gauge, d^A"^ = 0, where all the non-local terms vanish and ([4]) simplifies to the already 

mentioned operator J d'kc (A^)^. Let us also quote here that, in a generic linear covariant gauge, the anomalous 
dimension of A^^^ has been calculated at one loop order in [5T] . 



More recently, we have pointed out [52^ that another kind of gauge invariant nonlocal operator might be rele- 
vant, namely 

= Tr J d^xF^^iD^y'p^, . (5) 

Unlike expression ([3]) , the operator ([5|) has the advantage of being localizable by means of the introduction of a 
suitable set of auxiliary fields. Firstly introduced in the case of 3d YM 54J, the operator O has received renewed 
interest in the context of 4d YM. In fact, we have been able to show that, when cast in local form, it gives rise to 
a local action which can be proven to be renormalizable to all orders in the general class of the linear covariant 
gauges [SO]. Moreover, in the anomalous dimension of ([5]) has been evaluated at the two loop order in the MS 
scheme and explicitly proven to be independent from the gauge parameter. Also, in the case of the Landau gauge, 
it has been shown [53j that the inclusion of Zwanziger's horizon function does not spoil the renormalizability of 
©■ 

Despite the progress already achieved in the Landau and covariant linear gauges, a detailed analysis of the gauge 
invariant operator O in the MAG is still lacking. The main task of the present paper is to fill this gap, i.e. to 
achieve a local and renormalizable framework for O in the MAG. In particular, the manifest gauge invariance of 
the operator O might be useful in order to improve our present understanding of issues like the Abelian dominance 
and the dynamical gluon mass generation in the MAG. 

This paper is organized as follows. In section 2, by means of the introduction of auxiliary fields, we obtain a 
local action for YM theories in the MAG, in the presence of the mass operator ([5|). Moreover, the embedding 
of the resulting local theory in a more general action, enables us to make use of the BRST transformations. In 
section 3 we obtain the full set of Ward identities fulfilled by the starting action. Section 4 is devoted to the proof 
of the renormalizability of this action to all orders of perturbation theory. We obtain the most general invariant 
counterterm and we prove that it can be reabsorbed by means of a redefinition of fields and parameters of the 
starting action. The last section collects our conclusions. Appendix A contains a detailed discussion of the mass 
operator ([5]) in the presence of the horizon function for the MAG. 



2 Local action in the maximal Abelian gauge 

As is well known, the action ([T]) is left invariant by the gauge transformations 

6^A^ = -D^^u;^ , (6) 

for arbitrary u!^{x). In order to quantize the theory, one must fix the gauge. As we shall choose an Abelian gauge, 
we decompose the field as 

^ A^T^ ^ A;t^ + A^^T\ (7) 



3 



with 

[T°',T^] = igP^'^T" + igr'''T\ (8) 

[T^T'^j = igf"''T\ (9) 

[T\T^] = 0. (10) 

The index z = 1, . . . , iV — 1 labels the — 1 diagonal generators {T*} of the Cartan subalgebra of SU{N), while 
the index a — I, ... , N{N — 1) labels the remaining off-diagonal generators {T"^}. 

Accordingly, the field strength decomposes as 

f;, = D';!'Al~DfA''^+gr''^AlA:, (11) 

= d^Al-d^A^^ + gr'^A^Al, (12) 

where we have introduced the covariant derivative 1?°'' with respect to the diagonal components of the gauge 
field, namely 

= S-'d^ - gr'^A^ . (13) 



2.1 Gauge fixing 

The MAG [311 ISSl IMj is obtained by requiring that the functional 

n[A] = j d^x A^AI (14) 

is stationary with respect to the gauge transformations. Observe that expression (IT4l) depends only on the off- 
diagonal components of the gauge field. The vanishing of the first variation of TZ leads to the non-linear condition 

<4 = ■ (15) 

Still, it remains to choose a gauge condition for the diagonal components A^^ of the gauge field. We shall impose a 
Landau type condition, also employed in lattice numerical simulations, i.e. 

d^Al = . (16) 

The conditions (|15ll6p are implemented by adding the following gauge fixing term to the Yang- Mills action 



5mag = j d^x{b'^DfAl + c'^DfDl^c'' + gr''^c'^{Dl^A'^^e + 7fDf{gf^''Alc'') 

where 6"^ = {b°',b^) are the Nakanishi-Lautrup fields, and = {c°',c'), = (c",?) are the ghost and antighost 
fields. The gauge parameter a is introduced in (jl7p for renormalizability purposes. As a consequence of the non- 
linearity of condition psp . the quartic interaction ghost terms in expression ()17p is required in order to obtain a 
stable action [37l [56l [57] . After the removal of the ultraviolet divergences, the limit a ^ has to be considered in 
order to achieve (fTSl). 



In non-Abelian gauge theories one has to face the existence of the Gribov ambiguities [3], which deeply affect 
the infrared region. In the MAG, it is known that condition ([TSt does not uniquely fix the gauge [40], so that a 
suitable restriction of the domain of integration in the Feynman path integral has to be implemented in order to 
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avoid the counting of equivalent field configurations. The renormalization of the operator (O has already been 
investigated in the Landau gauge [53] when the restriction of the domain of integration to the so called Gribov 
region was taken into account. In [53j . it was explicitly shown that the introduction of the Zwanziger horizon 
function [11[5], which implements the restriction to the Gribov region, does not spoil the renormalizability of ([5]). 
The same feature occurs here in the MAG. However, for simplicity, we have decided of not including in the main 
text the lengthy and technical analysis of the mass operator ^ in the presence of the horizon function of the 
MAG, leaving the inclusion of the horizon term to Appendix A, where the interested reader may find a detailed 
discussion. 



2.2 Localizing the mass operator 

In [SD|, it has been shown that the nonlocal mass term ^ may be written in a local form by means of the 
introduction of a pair of complex bosonic antisymmetric tensor fields in the adjoint representation, (i?^^, B^^), and 
a pair of anticommuting antisymmetric complex tensor fields, {G^^, G^^). Following [50], the nonlocal operator O, 
eq.®, is coupled to the Yang-Mills action by introducing the gauge invariant mass term 

So = -^J d'x f;^^[{d')-']^'' f;^^ . (18) 

Furthermore, it is easily checked that expression (|18p can be rewritten in local form as 



e 



with 

Sbg 



j DBDBDGDGexp[-{SBG + Sm)] , (19) 
\ J d'x{B^^O^^B^^ - G^O^^Gi:) , (20) 



= ^d\{B%-B%)F^^^ (21) 

qAb ^ D^CjjCB (22) 

The identity ([TO]) allows us to localize the expression So, (|18p . when added to the YM action ^ym, eq.([T|). Thus, 
for the local gauge-fixed action S'phys in the MAG we write 

•Sphys — Sym + Smag + Sbg + Sm ■ (23) 

Evidently 

DBDBDGDGexp [-S'phys] = e-(^™+^"*<^+^°) . (24) 



and 



2.3 BRST invariance 



In order to establish the BRST invariance of the local action S'phys, we shall employ here the same procedure of 
|5l [52l [50l [38] , and we shall embed the action S'phys into a more general one. Following [5l [52l [50l [38] , we introduce 
the system of external sources 

and replace the term S.,„ in (|23p by 

Suv ^ 4 y '^'^^FpAUxppyGxp + VxppyBxp - VxppyB^^ + UxppyGxp) ■ (26) 
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Notice that expression (|2ip is recovered from (pS)) in the physical hmit for the sources, namely 

im 



U, 



I.e. 



Thus, action (|23p is replaced by 



phys 
phys 



= u, 



phys 2 

= 0, 



phys 



5i 



phys 



•Sym + 5'mag + Sbg + 'S't/v , 



which defines a more general theory and has S'phys as a particular case. 

As an important bonus, we observe that S'inv is invariant under the following global C/(6) symmetry 

Q/ii^Ap ('S'inv) — , 

with 



(27) 
(28) 
(29) 

(30) 



Q>pv\p 



Xp 



puctl3 



Uxpap - 



V, 



The symmetry (1301) is naturally related to the mass operator and allows us to use a multi- index notation, {pL,v) — > /, 
/ = 1, . . . , 6, so that 



V\pa(3 - 



(31) 



{Bf^Bf^Gj^Gf) 



-(B^ (^A pA \ 

X^f.LU- -^-^pw ^ pvi ^ pv ) 

I 



{Ulpv 1 U I , ^1 pu , ^Ipi/^ 2 (J-^crppi^ ' Uappu i ^(7pp,v i ^appu ) • 



(32) 
(33) 



The use of the multi-index / will turn out to be very useful when looking for combinations of possible counterterms 
respecting (15(11) . With this notation, we may rewrite the action (1^5)) as 



YM 



'5'mAG 



d^xiBfO^^ Bf - Gf O^^Gf 



Pau{Uj,_l_yGf + Vj,_,yBf — Vl^yBf + Ul^yGf)) 



(34) 

The introduction of the system of sources (pS)) enables us to establish that the action p4p possesses a BRST 
invariance. In fact, it can be checked by direct inspection that the following nilpotent transformations 



sc" = 

sc^ = 

sBf = 

sBf = 

sGf = 

sGf = 

sVif_,^ = 

sUipu = 



-idpC^+gr^'A^^c"), 

2 

6^ , = , 

gjABC^B^C 

Bf 



gfABC^BQC 

Uipu , sUi^p = , 
Vipu , sVfp^ = . 



leave the action p4p invariant. In particular, we can rewrite 
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MAG 



{ Df 



abi -a -b i 



rabc a-b -c 

-t9J c c c 



(35) 
(36) 
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3 Identification of the complete starting action and its symmetries 



In order to prove the renormalizability of the action 5i„v, ea. (|34p . we shall make use of the algebraic renormalization 
technique [58]. To that purpose, we shall introduce a set of BRST sources and we shall establish the Ward identities 
fulfilled by the theory. The knowledge of these Ward identities will play a central role in the determination of the 
most general allowed invariant counterterm. 



3.1 Starting action 

As is well known, due to the nonlinear character of the BRST transformations, eqs. (|55|) . one has to introduce a set 
of external sources, {Xf, Xf, Yf^ , Y/^, Q^, L^), coupled to them, namely 



(37) 



+ ^pU^^^a^b ^ Yj^gjABC^B^C _^ y/" g f^^^ c"" Bf + Xf g f ^^^^ G'i + Xf g f ^^^^ O'i 
Further, we introduce the following two terms, S\ and S'souiccs to be added to the action Sim'. 
Sx ^ J d^x S^Xi{BfBf - GfGf){Vj^,Vj^, - Uj^^Uj^,) 

+A2 (^fG^Vi^uUjfji,v + GfC^Ui^uUjiiu + BfB^Vi^uVj^u — GfBfUi^,yVj^i, — CfB^Ui^^Vj^j^^ 
+GfBfUi^,jVj^ij — -BfB"fVi^uVj^y + -GfG'jUi^uUjfj.i, — -BfBfVif_,„Vj^i, + -GfC^Ui^yUj^u 



\ABCD 

-{BfBf-GfGf){BSB?-G3G'^)\, 
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d^x 



C{Ulf_,i,Ulf^,yUjapUjap + Vl^vVl^vVjapVjcp — "iU I p^U I pvVj a pVj a p) ■ 



(38) 



(39) 



The first term, already discussed in detail in |50[ 155] . contains interactions between the auxiliary fields and the 
sources {Uipv,Uip^,Vip^Vip^), and is needed for the stability of the action. The second term, S'souiccs, depends 
only from {Ujp^, Ujpi,, Vipi,Vipi,) and is allowed by power counting. The parameters Ai, A2, Xi, X2 and C are free, 
while the 4-rank invariant tensor \^bcd enjoys the properties [SOI 



MAN \M BCD I fMBN \AMCD , rMCN \ABMD , fMDN \ABCM 



J MAN ^Jl 



jMiJN -^ABMLI _j_ JMJJN 



= 



and 



A 



ABCD _ yCDAB _ ^BACD 



Collecting all terms, we finally write the starting action S as 

^ — Sinv ^" S\ 4" Scxt ^~ ^sources ■ 

The stability of S under quantum corrections will be investigated in the next section. 

Here, for the sake of clarity, we recall the range of variations of all the indices introduced so far; 



(40) 
(41) 

(42) 



A,B,C,D,.. 


e 


{!,• 


.,iV2-l}, 


a, 6, c,d, . . 


e 


{!,• 


.,iV(7V-l)}, 


h jj k^l, . . 


e 


{!,- 


.,iV-l} , 


I,J,K,L.. 


e 


{!,- 


.,^D{D -1)^6} 


(J,p,.. 


e 


{!,• 


.,D = A} . 



(43) 
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A 


6 


c 


c 


B 


B 


G 


G 


[/ 


U 


1/ 


V 


Y 


Y 


X 


X 


n 


L 


dim 


1 


2 


2 





1 


1 


1 


1 


1 


1 


1 


1 


3 


3 


3 


3 


3 


4 


gh. number 








-1 


1 








1 


-1 


1 


-1 








-1 


-1 


-2 





-1 


-2 


Qe-charge 














1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 









Table 1: Quantum numbers of the fields and sources 

In table ([T]) we display the dimensions and the ghost number of the complete set of fields and sources of the theory, 
as well as the corresponding Qg-charge, which is defined as the trace of the of the operator ([31]) . 



The starting action (|42|) is left invariant by the action of the nilpotent BRST operator s given by ([35]) and by 

sO;^ = , sL^ = , 

sy/ = Xf , sXf = , 

sXf = -Yj-^ , sF/ , (44) 

i.e. 

sS = . (45) 

3.2 Ward identities 

In order to constrain the possible counterterms which can be added to the local starting action S, eg . (|42| . let us 
proceed by establishing the set of Ward identities fulfilled by the starting classical action. In fact, it turns out that 
S obeys the following set of Ward identities: 



• The diagonal gauge fixing identity: 

• The diagonal anti-ghost equation: 

• The diagonal ghost equation: 
where 



§ = C.A;, . (46) 



g'iE) = A^i_ , (48) 
a^ = ^+./-c-|,, (49) 



and 



A^iass = -d^c^~d^n'+gr>''niA''^gr''^L-c''+gr''%-B'j 



MM 

Notice that the term A^j^^^g, being linear in the quantum fields, represents a classical breaking not affected 
by quantum corrections |58| . 

• The Slavnov- Taylor identity 

cS(E) = , (51) 
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where 
5(E) 



f 5Y. 



(5E (51] 



5Y. 5Y. 



Bf 



(5E - (51] 



(5E 
(5E 



Gf 



(5E (5E 



(5E (5E 



(55^ (5y/ (55/ (5X/ SGf 



(5E 



(5S 



The Slavnov- Taylor identity (|5ip gives rise to the corresponding nilpotent hnearized operator 



( (5E 



(5E 



(5E 



(5E ^ (5E J 



i^^&Af; (5A;4^f^A J^A^^A ^c'4(5L'4 



V(5F/ 
(5E (5 



(5E J 



(5E 



(55/ (55/ (5r/ (5F/ (55/ 



(5E 



(5E 



SXfSGf SXfSGf \SGf 
'SGf SXf^ ''"'SU, 



(5E 

Zb/ 

(5 / (5E 

M/ ^ vZx/ 



(5c^ 



5/ 



• The diagonal C/(l)^-i Ward identity: 
with 



5s5e = . 



SA] 



+ ^^r^ ( ^r.^ + + ^'^^ + + 5?^ + 5?^ 



(55| 



' 5G] ' 5G\ 



5_ 
)B 

5 



SB'} 



(5 

SGf 



XI 



sxi 



(52) 



(53) 
(54) 
(55) 



(56) 



• The off-diagonal SL{2,R) identity: 
where 



I?(E) 



P(E) - , 

d X C 



Sc" 5L°- 61 

As the Slavnov- Taylor identity, the equation ([57]) also defines a linear operator: 

JE (5 



j4 / a ^ (5E (5 

d^a; c°- 1- Tl—Tl- + Tl-TT- 

Sc'' 6h°- (56" (5L" 



(57) 
(58) 

(59) 



• The global C/(6) invariance related to the nonlocal mass operator: 

Q/,/(E) = , 



where 



Qij ^ d'^xl Bf 



6Bf 



SBf 
5 



^'5Gj'^-' 



SGf 



~Y 



Xi 



(5 



x\ 



'5Vjt,, ""^"(5^/^. ■ ' 5Yf (5y/ ■ 5Xf 5Xf 

The trace of (|6T|) defines the Qe-charge, Qq — Qu, already displayed in table H]). 



(60) 



(61) 
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nfj\i:) = Q, (62) 



• The exact rigid symmetries: 
being TV = 1, 2, 3, 4 and 

nf] EE jd^.{Bf 

T^u ^ J d'x {sf^ + Gjj^^ ^/.^TTT— ^ Uj,,— - ^/t^ + ^fj^ ) ■ (63) 

For calculation purposes, we display the following useful (anti-)commutation relations between the linearized op- 
erator Ss, given in ([55)1 . and the operators (|^ . and ([55]) . namely 

2?s,5s]=0, {g\S^}^W\ (64) 
[n\'j , 5s } = Q/j , {Te^'j , 5s } = , (65) 

{7^fj,5s} = 1 A (<5,;,<5^i - <5zi<5,7K) (^^^ " ^^^^^ + ^) ' ^^^^ 



6 


-Gj 


S 


+ Vif,^ 


s 




s 


+ r/ 


s 


^xj 


<5 


SGj 


SBj 


SUj^^ 






sxj 






S 


+ Gj 


S 


+ Vif,^ 


s 


+ Ujf,^ 


s 


-y/ 


s 


^xj 


5 


SGj 


SBf 


SUj^u 


SVi^,, 


sxj 




SYf' 


s 


-Gj 


S 




s 




5 


+ F/ 


s 


^xj 


5 


SGj 


SBf 


SUj^u 


SVi^, 




sxj 




SYf' 


S 


+ Gj 


S 




s 




5 




s 




5 


SGj 


SBf 


SUj^^ 


SVi,,, 


sxj 




SYf' 



is: s: , , , \ i r^A 5 jj ^ ^ 



d'x (SikSjl + SilSjk) [G'kj^ - f^^A.- '^^^) ■ ^^^^ 

4 Renormalization 

In the previous section we established the full set of Ward identities fulfilled by the action S. In this section, we 
prove that E is perturbative renormalizable to all orders. 

4.1 Determination of the most general counterterm 

Let us turn our attention to the characterization of the most general invariant counterterm Sqt which can be 
freely added to S. According to the algebraic renormalization procedure [58] . we require that the perturbed action 
(S + sEct) satisfy the same set of Ward identities, and constraints of E. The counterterm Ect must be 

an integrated local polynomial in the fields and sources with dimension bounded by four, vanishing ghost number 
and Qe-charge, obeying the following constraints 

5e(Ect) = 0, (68) 

JEcT 



, (69) 



2?e(Ect) = 0, (71) 

QK^ct) = 0, (72) 

W^(Ect) = 0, (73) 

Q/j(Ect) = 0, (74) 

n\^j\^CT) = 0. (75) 
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As an immediate consequence of BRST invariance, condition (|68p allows us to write 

ScT = ao Sym + 5eA(-i) , (76) 

where Oq is a free parameter and A*-"^) is an integrated local polynomial with ghost number —1 and vanishing 
Qg-charge. Taking table ([1]) into account and imposing the conditions and (|70p. we are led to the expression 



with 



A(-i) = A^-^) + a'-^) , (77) 

A*-') ^ J {«i + a2 (9,.c")A;J + ag {^], + a^c*)^ + 04 c"L" + 05 c'L* + ae gP^'^cfc^e 

+ar gr''^c''c''c' + as c'^b'^ + ag gr^'c'^A'^Al + Y/^Bf + an Y/^Bf + ai2 Xfcf 
+ai3 xfGf + ai4 Gfd^Bf + ais gf^''''Gf{d,A^)Bf + a,e gf^'^'^GfAf^d^Bf 
+ai7 {d^A^)Ui^,,Bf + ai8 Api^^.d^Bf + a^g gf^""^ A^A^Ui^,,Bf + a^o {d^,A^)Vi^,Gf 
+021 A^Vi^,,d,Gf + 022 g/^^^A^Al; V^/^.Gf + a23 Of BfVj ^,,Vj + 024 QfBjVi^.Vj^, 
+025 GfBjVj^i,Vi^i, + a26 Gf Bfu J ,_i_yU J + 027 GfB^Uj,_i_yUj,_i_y + 028 Gf Bfu j ^,^11 i 
+a29GfGjVi^i,Ujf_ii, + a3oGfGjVji_i,^Uif_ii, + a^nBf B'^Vi ^,yU j + a32Bf BfVjf^^Ui^i, 
+033 BfBfjJjf^uVj^u + 034 BfBjUif^^Vjfj_i, + 035 GfGjUifj^^Vjf^u + 035 GfBfVjf_,„Vjfj_i, 
+037 G^^^:^^/^.K7m. + ^38 GfBjUi^Mj^.. + P^^'''' GfBfA'lA^ + 039 c'^c'^C/7M^'^/M^' 

+^fcd-a^bQc^d ^ ^aMj ^a^bQ^^^j^ _^ ^ABCD qA _^ ^ABCD qA^BqCqD 

+a3^^^^ GfB!^BfB^ + (74^^^^ G^^i?f Gf G^*} , (78) 



and 



A2 



(-1) _ 



d'^X ^040 C//pj/<9^V7pi/ + 041 Ulf^^df^daVl^a + 042 Ul^^Vlf_i^Ujaf}Uja/3 

+043 ?7/^i,V7^,^VjQ^Vja/3 + 044 Ul^^V,Jfj,^Ula0U,Jap + 045 V^/^i/ Vfa/jVja^ 

+046 Uj^^VluaUjapUjfi^i, + 047 CZ/^tiy Vfi^Q K/q^ + 043 U I ^pVj yaU lafjU J 

+049 i7/^j^Vji,aVfa/3Vj/3^ + 050 Vf^j. Vf^i/J/ja^K/a^ + 051 Vl ^yVj ^yU JafiVj aji 

+«52 Ul^yUj^yVjajiUjaf} + 053 U I ^^VlvaU J ,_,pU J af} + 054 U I ^vVlvoiVj ^ijiVj afS 

+ a55 Ul^yVjvaUl^pUjaP + 056 CZ/piyK^i^a^/Ai/sK/a/S + 057 U J ^i,VlapU J apU J fiy 

+ 058 C^/^i/V/q/jFjq/jVj^^ + 059 UiafiVlvpUjal3Ujf_,iy + OgO t//aM^/i'/3K7Q/3K7/jiy^ (79) 

Here 0^, a = 1, ■ • • , 60, and the four-rank tensor a^^'^'^ , K = 1, . . . , 4, are also arbitrary coefficients. Applying the 
remaining conditions ([7TH75)l . expression (|77|) reduces to 

A(-i' = y d^xjai - as c'^L'^^AI^ + 04 c"L° + ag (c"6" - gP^'c^c^d - ^f'^^c^c^c"^ 

+aio (r/B/ - X/Gf ) + an {Y.^Bf Xfcf) + au CfofO^^^Bf 

+019 F^yUi^ijBf + (flio — On - flig) F^j,V/^,^Gj + 023 GfBf (K/^^Vj^^ — Oj^j^yUj^^) 

+025 {GfBjVj^i,Vi^i, + GfGjVif^i,Uj^i,) + 028 {Cf Bjtj J ^ijU I + Bf B^Vi ^lylJ J 

+ |[A2(aio — On) — 028 — 025] {Bf BfU i ^^Vj ^j^y + GfBfUi^yUj^y) 

-^[A2(aio - On) + 028 + 025] {GfGjUi^,yVjf,y + GfBjVif,yVj^,y) 

^^ABCD QA^B ^^C^D _ ^C^D) ^ CZ/^.aV/^. + 041 ^J/^.^M^a^/.a 

+O42 (^//^i^F/^j/C/jQ/jC/ja^ — ?//^i/V7^^Vjct/3Vja/3) I , (80) 

with the further restrictions on af^'-^^ 

rMAN MBCD I rMBN AMCD , rMCN ABMD , rMDN ABCM n /Ol^ 
i CTl + / cr^ + / (7^ + / (Ti = , (81) 
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and 

ABCD CDAB RACD fQc,\ 

Observe that many coefficients present in expressions ([751) and (|7^ vanisti after tlie imposition of tlic constraints 

(tZTHZS- 

Performing now the foUowing redefinition 

a\ — > ai , 

02 -02 , 

04 -03 , 

a 

fflio + an as , 

aiO — 019 ^ 16 : 

023 ^ A1O7 , 

025 + 028 ^ A20g , 

040 X1O9 , 

041 ^ X2O10 , 



042 Can 

ABCD a^(„ABCD \ABCD 

'1 
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we can rewrite the counterterm Ect as 

ScT = EcT(l) + EcT(2) + EcT(3) (84) 



12 



ao + 



' ao + 2ai 



(DfAi) {Dl^Al - Dl-A'i) + (ao + M^jgr"" {DfAt) A^Al 



4ai 



/ 

~l~ ^ V-^ •! " "/i" "t'"' "/J." "Z^ ' -"^ " '/.i" "/V '-fl" '"h' J ' 2 

+ (ao + 201)5/""* (a^^t) A;a1 + {a, + a^^DfA; + {a, + a^)gr'^^lf (i^^^A^) 
+ (ai + 02 + 03)0^^7^^ {gf^^Alc") + (2ai + 02 + 03)3 VV^c'^C^^^^ 



-(ai + 03) (1^; + a^?) gr'M^lc" + (02 + a^)<fDfDl-c^ + a^^b'^b'' - a^agr'^'b^'^d 

-(«3 + ^4) [^gr'^b'^c'^c^ + J^'r'v" 



,2 ^abc ^ade^^c^d^e 



C C C C 



(2a3 + 04) 



+ (ai - as) niDfc' - as gr'^i^^A^c' + as |/»"«L»c"c^ 



+as5r'''i'c"t 

y d^x I as {gr'"c''B'j + gr'^c'B}) + as r^^Y^c^B^j + as F/ {gr'"c''Bj 

+ gr^'c^B}) + as gr^'Yjd'B] + as X," {gr^'c''Gj + gr^'c^G}) + as gt'^Xic'^G] 

+as {gr'"'c'G'i + gr^'c^G]) + as gr'^Xlc'^&'j + a, B'^d^B'} - ag 5/"''^? (9^4) 

-(«i + a5)gr'"Bf (a^4) - as 2gr''B]A'^d^B'j - (ai + a5)25/'"'^5?A^a^Sj 

+a5 5'r''r''''5?44^? + («i + a^)gT'r'"'B<iAlAlB'l + (ai + a^jg^r^T^'B^Af^A'^B 

+(2ai + a^)g^f'""^r^''B'iA'lAlB'i + (2oi + a^)g^ T""' f'^' B'} A^A'Ib] + 05 Sja^B} 

-(2ai + a5).9^/"'\/""^S}A^4i3| + (ai + a,)gr''B] (9^4) + (ai + a,)2g r"' B^ A%B 

-{ai + a^)g'^r'^ r'^'B'iA^^AlB] - (2ai + a^)g'^ f"^' B^j A^A^B] 

-{ai + a5)gr''& (9^4) Bf - (oi + a^)2gr'>'B\Ald^B'i + (ai + a^)g^ r'T' B^A^A^B'^ 
+(2ai + a5)5'r'7'='"5}44S?| , 
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and 



^CT(3) 



-(ai + a5).g/'"'^G? (a^A^) GJ - as 2g/'^''^G?A;9^Gj - (ai + a5)2g/''''=G?A^a^G5 
+a5 5V""/^''^G?A;A;,G| + (ai + a5)5V'^^7'''G?A^A^G| + (ai + a^)g^ r^'G'iAlA'^G] 
+ (2ai + a5)gV°'''/''"G?A;:A^Gj + (2ai + a^)g^ f'^^G'^AI^AlG] + ag G^a^G} 
-(2ai + a5)gV°'7°'^G}4A^Gj + (ai + a5).g/"''^G? (9^4) G} + (ai ' - ^o„.ab»/.a 
+ a^)g\n r'^G'iA^^AlG] - (2ai + a,)g^ T" f'^G'iAlA'lG] 



a,)2gr'^G^jAid^G] 



+ (2ai + a5)gyT''G}AlA'lG'}\ + [2(ai + ag) i?;^^^^ + (2ai + ae)^/"''^^^^] (C7,^, G? 

+ V7^. - V>M. + C//^, G?) + [2a6 + (2ai + ae)gr''A;Al] (Uj^, G] 

+ Vi^,, B] - Vi^,, B] + Ui^, G\) + Ai(a5 + ay) {BfBf - GfGf) {Vj^.Vj^, ~ IJj^.Uj^,) 

+A2(a5 + a%){j3fBjVi,_i_yVj^y — GfBjUi,_i_yVjf^i, — Gf BfUif_i^Vj^i, 

+BfGjVif^„Ujf_ii, + GfGjUif^ijUj^u + GfBfUifj_,^V.j^„ — ^BfBfVifj_^Vjfj_i, + \GfG^Ui^ytJj^y 



(«i + ci,)gr'^& (9^4) G? - (ai + a5)25/"^'Gl49^G? 



ah Ah 



(ai 



a^)g^r'^r^G\AlA=^G'i 



^itiuVjfiu + \Gi GjUi^i,Ujf^,y ) + —(A 



"5 /\ABCD 



ABCD 



){BfB 



I -GiGi){BjBj -urjLTj 



Cr<D\ 



~C,aii{Ul^uUl^vUja(3Ujal3 + Vj ^^Vj ^^Vj apVj ap " "^Vl ^vVj ^^U J ajiU J ap) 



G\G 



(87) 



By construction, expression 
identities. 



yields the most general invariant counterterm compatible with the full set of Ward 



4.2 Renormalization factors 

Once we have found the most general counterterm, eq. (|84p . we have to check whether the remaining independent 
coefficients, oq, ai, . . . , an, and the 4-rank tensor cr-^sc'D reabsorbed through a redefinition of the fields, 

sources and parameters of the starting action S. The answer is in fact affirmative. Let us rename collectively the 
fields, sources and parameters as: 



$ = (A, 6, c, c) and ^ = {B, B, G, G) , 
J={n,L) and J ^ {Y,Y,X,X,U,U,V,V) , (89) 

f (5,a, Ai,A2,xi,X2,C) • (90) 
Next, defining the bare fields, sources and parameters as: 

^oflf-diag _ yl/2 ^off-diag 
J off-diag _ Z J J of^'^ii^g 

<Jo = Zj J , 

Co = ^eC, (91) 



and 



yABvu ^ A^^^^ + Z"""^^ , (92) 
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it is easily checked that the invariant countertcrm Sct can be reabsorbed into the starting classical action S, 
namely 

M/, J, J, + e EcT = SKo, *o, *o, Jo, Jo, A^^^^] + 0{e^) , (93) 

where e stands for an infinitesimal expansion parameter. For the Z's factors we have 



with, 



A 


= z-^ , 






= ZaZ^'^Z^'^ 




zV 


— z 




z-J 


= zy^ , 






_ 7-1/2 




yl/2 


- - 


= z, 


Zn 


= z;'z-'/'z- 


1/2 

) 


Za 


= z-^i"^ , 




Zl 


= z-'z-^i^~z- 


1/2 

) 


Zl 






Zx 


= Zx = Zy — 


Zy = 


Zu 


= Zjj = Zy = 


Zy , 



1/2 
B 



7-1 7-1/2 7-1/2 



71/2 








Zg 








zV" 




, , 02 + 03 

' + ^ 2 ' 




zT 








zT 




1 + ey , 




Zy 




1 /Oo 05 




Za 




1 + e (oo - 202 - 


1-04) , 


Zx, 




1 + e (oo + 05 - 


206 + 07) , 


Zx. 




1 + e (oo + 05 - 


2a6 + as) , 


Zxi 




1 + e (oo + 05 - 


206 + ag) , 


Zxi 




1 + e (00 + 05 - 


206 + 010) , 


Zc 




1 + e (2ao + 2a5 


- 4a6 +aii) , 


Zx 




1 - e 05 , 




■ABCD 




ABCD 





(94) 



(95) 



This concludes the proof of the renormalizability of the classical action to all orders of perturbation theory. 



5 Conclusions 



In this paper, a detailed analysis of the nonlocal gauge invariant mass operator Tr J d'^xFij,i,{D'^)~^ F^^ has been 
made in the MAG. By means of the introduction of a suitable set of auxiliary fields, this operator can be cast 
in local form. Moreover, the embedding of the resulting local model into a more general action has allowed us 
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to make use of the BRST symmetry. Furthermore, it turns out that the generahzed action displays additional 
global symmetries giving rise to useful Ward identities, which were used to restrict the possible counterterms. The 
analysis of the renormalization factors has enabled us to show that the most general invariant counterterm can 
be in fact reabsorbed into the starting action through a redefinition of fields, parameter and sources, establishing 
thus the perturbative renormalizability of the model to all orders. Finally, in appendix A the nonlocal operator 
Tr j d'^xF^^[D^)~^ F^i, has been analyzed in the presence of the horizon function implementing the restriction of 
the domain of integration in the Feynman path integral to the Gribov region in the MAG. The output of our 
analysis is that the introduction of the horizon function does not spoil the renormalizability of the model. 
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A Including the horizon function 

It is a well known fact that non-Abelian theories are plagued by Gribov ambiguities ^ , see f60^ for a pedagogical 
review. In the specific case of the MAG, the study of Gribov ambiguities and the characterization of the horizon 
function in the case of SU{2) can be found in [39l[38]. In this appendix we present the main aspects of the simul- 
taneous inclusion of the MAG horizon function and of the gauge invariant mass operator ^ in the YM theory. A 
similar treatment regarding Gribov ambiguities and the mass operator ([5]) has been done recently in the Landau 
gauge |53J. Without loss of generality, we will follow [39l[38] and restrict ourselves to SU{2). 

In [3S] , the following horizon function for the MAG has been derived 

^Horizon = 1^9^ J d\ E^" . (96) 

Here 7 stands for the Gribov parameter [31 = e'^"'', a,b ~ 1,2 are the off-diagonal components of the 

SU{2) structure constants and A^^ = A^ is the diagonal component of the gauge field. The operator (A^"^)"'' is 
the inverse of the Faddeev-Popov operator given by 

^ab ^ -Dl'^Df - g'^e'^'e^'^AlA'l^ , (97) 
with the covariant derivative defined as a particular case of by 

Df = S-^'df, - ge'^'A^ . (98) 

The inclusion of the horizon function (|96p allows one to implement the restriction of the domain of integration in 
the Feynman path integral to the Gribov region, where the operator (|97p is strictly positive definite. As underlined 
in [321 [35], such a restriction is necessary in order to deal with the Gribov copies. 

In much the same way as the gauge invariant mass operator ([S|) , the horizon function (|96p also possesses a localized 
version, which can be obtained through the introduction of a suitable pair of commuting auxiliary complex fields 
{<t>'^,<t>'^), and a pair of anti-commuting ones (w°'',a;°^) [3H] 



(99) 
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with T given by 



(100) 



and the sources Af!J^, M^^, A^°^, and Nf^^ are chosen in such a way that their physical values must be taken as 



m; 



prab 



phys 



phys 



= -M, 



cab c z 
= Of_i^J , 



phys 

= 



phys 



(101) 



Notice that the last term in expression (I99p introduces quartic interactions between the auxiliary fields, being needed 
for renormalizability. Nevertheless, unlike the term psp . these quartic terms depend on the gauge parameter a of 
(fT7|) . which will be set to zero, a ^ 0, after the removal of the ultraviolet divergences. 



Adding the local version of the horizon function to the action (|^^ . we obtain a new starting action which 
reads 

"5* — Sinv + '5'hortzon ^" '^'sources + "5*0x1 • (102) 

The first three terms of expression (|102p are given by and ([55)1 respectively, while the fourth term 

generalizes including the new sources (|10ip . namely 



-•sources 



5. 



sources 



j^abpjab 



(103) 



Finally, the last term of (|102p contains the couphng of the external sources fij^, r°, ^JJ, fi^, L", L, ,YP , Xf ,Xf 
and , rilf , pif , 'dif to some nonlinear operators needed for the BRST invariance of the model, being given by 



+ ge^'pl^Col^c-dl 



ge^b^be^ _ 9_^ab^cd-b^e^c^d 



(104) 



The quantum numbers of the new fields and sources are displayed in table [21 

Let us now proceed by giving the set of Ward identities fulfilled by the action (|102p . These are: 

• The Slavnov- Taylor identity: 



SiS) 



SS 5S \ 6S 5S 6S 



Sn- St- J 6 A- Sn^SA^ 



6S SS 5S 
SlAS^^ sE^ 



SS SS 



SXfSGf \SXf 



SS 



Bi 



SS - SS 



jab_^_^lab_^^^ab 



SGj 
SS 



SUi^^ 
SS 



SS 



SS 



Y 



Gf 



SS 



SS 



SS ST. 



5Bf SYf SBf 
A SS 



X 



M°''^ — 



''""svi^, sxf "^"^ <5r/ 

SS 



tna SS ^gfc SS 



S^f 



. (105) 



The global t/(8) invariance: 



with 



QtuiS) - , 



(106) 



Q 



b _ 



Icb 



i,cb 



'sm 



""SN: 



S 



Si}?!' 



Sn 



ca 



\ cb 



Spf 



(107) 
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u M 


M 


iV 


iV 


i9 77 


P 


A 


n 


r 


C 


dimension 


1 


1 


1 


1 2 


2 


2 


2 


3 3 


3 


3 


3 


3 


3 


gh. number 








1 


-1 





1 


-1 


-1 -2 





-1 


-1 


-1 


-2 


Qg-charge 


1 


-1 


1 


-1 1 


-1 


1 


-1 


1 -1 


1 


-1 












Table 2: Quantum numbers of the fields and sources 

The presence of this global invariance U{8) allows us to make use of a composite index i = {a,fj,), with 
i = 1, . . . , 8. Thus, from now on, we set 



[<Pti ) 



(W^ M"'' N"^ iV"'') 



The trace of (I107P defines a new Qg- charge whose nonvanishing values are displayed in tabled 
• Symmetries involving the Faddeev-Popov ghost fields and the localizing fields: 

SS 



,5S_ 
5S 



i9f 



(?5 

SS SS SS 
Sc'' (5?7f (56" 
(55 



(55* 


SS SS' 




(55 


^5 SS' 




(55 


~pA — 


(5^,f 


' ' ) SL'^ 



= 0, 
= , 



= 



(55 



Spf J Sc'^ 



'''srif, 



SS SS 
&dfSb^ 
SS 



Sujf J SL" 



= 



(108) 
(109) 
(110) 



(111) 
(112) 

(113) 

(114) 



• The exact rigid symmetries associated to the horizon function: 



SS 
5^ 



d'^x Luf 



SS 



SS 



to A 



SN°:. 



„ SS 



til 



Sujf 



5pt 
SS 



= 



(55 



i9? 



S^ 



(5A? 



SS , „ (55 



(5i?? 



0, 



(115) 
(116) 
(117) 



The global U{6) invariance: 



Qij{S) = 



with the operator Q/j given by ((6T|) . 
• The exact rigid symmetries associated to the mass operator: 



with the same operators Tlf^^ , A/" = 1, 2, 3, 4, already defined in ((63|) . 



(118) 



(119) 
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The diagonal U{1) Ward identity: 
with 



W^{S) = -d% , (120) 



h 7V°-^ 1- 77°— + d""—— 















h Af 



















(121) 



The off-diagonal SL(2,E) identity: 



The diagonal gauge fixing: 



P,5,= /A(e-§ + §§)=0. (122) 



^ = i>fK ■ (123) 



The diagonal anti-ghost equation: 



SS_ SS 



d.jT^^O. (124) 



We notice here that no diagonal ghost equation, similar to (pS)) . holds when the horizon function is taken into 
account [35] ■ However, the set of Ward identities listed above forbids the presence of counterterms like 

d'x {4>1^1 - Cjtu:t){BtBf - GfGf) , (125) 

as well as of any other counterterm which would mix fields associated with the two different nonlocal operators 
(O and ([TO)) . Therefore, in complete analogy with the case of the Landau gauge [S3], the two operators ([5]) and 
(j96[) do not mix, due to the rich symmetry content of the resulting local action. Moreover, it turns out that the 
most general allowed counterterm can be in fact reabsorbed in the starting action, (|102p . through a redefinition 
of fields, parameters and sources, ensuring the renormalizability of the mass operator ^ in the presence of the 
horizon function 
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